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ABSTRACT 

Recent Sternberg astronomical institute activities on site testing programs and technique are presented. The 
main attention is paid to the new modifications of MASS and DIMM data processing developed by SAI team. 
Four important unresolved questions affected to optical turbulence measurements are raised in the hope to be 
solved in nearest future. 
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1. INTRODUCTION 

In spite of that SAI site testing researches have been started many years ago, this presentation responds on 
our activity in last five years only. First, we present a short information about finished campaign for optical 
turbulence measurement at Mt. Maidanak and on-going site testing program at Mt. Shatdzhatmaz, where the 
new 2.5 m telescope is planned to be install. 

More detailed discussion is devoted to MASS and DIMM measurements processing, where some additional 
effects were taken in account in last years. Then, our plans for further enhancement of MASS and DIMM 
instrumentation and software are highlighted. Finally, some extra-essential questions are formulated in hope to 
be solved. 



Our team website is located at http://curl.sai.msu.ru/mass A lot of documents related to the MASS and 



DIMM instruments and corresponding software can be found on the site. We pay main attention to the both 
methods details. Also, SAI Automatic Seeing Monitor (ASM) is available at http://eagle.sai.msu.ru/ link. In 
the now presented works, Nicolai Shatsky, Olga Voziakova, Sergey Potanin, Boris Safonov, Matwey Kornilov 
took part. 

2. 2005 - 2007 CAMPAIGN AT MOUNT MAIDANAK 

Main goal of the campaign at mount Maidanak was to finalize 1998 - 1999 studieJi' of the optical turbulence 
(OT) above the summit. Our secondary intention was to extend our experience in MASS observations and 
data processing. The campaign was performed in collaboration with Tashkent astronomical institute staff. The 
observations have been started in August 2005 and finished in November 2007 have covered 5 seasons. 

Original MASS device installed at astrographic refractor was used as instrument. The telescope has aperture 
D=230 mm and focal length F=2300 mm. The original MAS^ was buift in 2002 in the cooperation ESO -f 
CTIO + SAI in 3 copies. The instrument is driven by Turbina software working under GNU/Linux. 

Data set was collected for 280 nights with whole length of 1022 hours. Data processing has resulted more 
50 000 OT vertical profiles. The main parameters of atmospheric OT are: free atmosphere (0.5 km and above) 
seeing Pfree — 0.47", isoplanatic angle 9o = 2.19". When seeing is better than its median then Oq = 2.47". The 
time constant is tq ~ 3.94 ms. Under weak turbulence tq — 5.41 ms. Corrected atmospheric time constant w 7 
ms. 

The results of this campaign are published in recent papeif^ where the prospect of adaptive optics system on 
the Maidanak 1.5 m telescope is discussed too. 
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Figure 1. Left: The map 2x2 km of the local relief. Right: Night time view of the ASM, July 2008. 

3. 2007 2010 CAMPAIGN AT MOUNT SHATDZHATMAZ AT NORTHERN 

CAUCASUS 

Mt. Shatdzatmaz (2127 m) is located in Karachay-Cherkess Republic of Russia, 20 km southward from Kislovodsk. 
The mountain belongs to the Skalistiy ('Rocky') ridge which is parallel to the Main Caucasus ridge «50 km away 
to North. Solar station of Main astronomical observatory (Pulkovo observatory) is situated about 1 km from 
the place chosen for installation of new 2.5 m telescope on the top of Mt. Shatdzatmaz. The nearest part of the 
local relief is shown in Fig. [T] on left. 

The main goal of this campaign is to collect statistically reliable data on seeing and OT vertical distribution. 
In parallel, the representative information on the amount and quality of clear night sky, on atmospheric trans- 
parency, sky brightness and on-site weather parameters had to be accumulated. This information is needed to 
develop the optimal strategy for the 2.5 m telescope operations. 

The ASM tower is installed at 40 m to SW from the spot reserved for the 2.5-m telescope. The ASM 
telescope tube is raised at 6-m elevation above the ground. Seeing monitor includes next instrumental and 
software components: 

• MASS/DIMM device to measure strength and vertical distribution of optical turbulence; 

• Telescope Meade RCX400 12" to feed MASS/DIMM; 

• CCD finder/guider to help accurate pointing to target stars; 

• Control computer to serve MASS and DIMM data acquisition; 

• Automated enclosure of close cloth; 

• Wind direction and speed, air temperature, air humidity sensors; 

• Boltwood clouds sensor on Solar station, self-made clouds sensor at the ASM dome; 

• Two web-cameras for internal and external overlook; 

• Controlled power supply for instruments; 

• Service computer to monitor environment and support ASM work; 

• Wi-Fi bridge to Solar station where ASM server is located; 

• Main software: mass, dimm, rcx^scope, monitor, dome, ameha; 



More detailed description of SAI ASM is available in recent paper. ^ Main properties of tested site is listed below. 



Figure 2. MASS/DIMM instrument installed on Meade telescope. The DIMM camera is Prosilica EC650 with IEEE-1394 
interface. Telescope finder/guider equipped with CCD camera is seen on the right side 

3.1 Clear skies and meteo-characteristics 

An estimation of clear skies has been done on the base of sky temperature measured by IR sensor: Boltwood clouds 
sensor. We have defined empirically that the sky can be "photometrically" clear if AT = Tsky — Tamb < — 22°C. 
For the site annual clear astronomical night skies w 1340'* or 46%. The maximum of the clear skies amount is 
observed from mid-September to mid-March, where about 70% of the clear weather is concentrated. 

Median temperature over year is -fl.8°C, temperature span is not very large: from +15°C in summer to 
— 15°C -15 in winter. Median wind speed is 2.3 m/s, the dominant winds direction is from west or south-east. 

3.2 Statistics of the campaign 

The total duration of the observations > 2700'* in the period 2007 November - 2010 August. Number of 
accumulated profiles > 130000. During the period the telescope has produced « 3300 pointings to target stars. 

ASM efficiency (used time to clear sky) « 50% in total. Over last year the efficiency > 75%. In 2009 two sub- 
programs were added: photometric program for atmospheric extinction determination and twilight observations 
of the OT. 

3.3 Main characteristics of the site 

At Mt. Shatdzatmaz overall median seeing (3 = 0.93". The most probable seeing value is 0.82". In 25% of the 
time seeing is better than 0.73". The free atmosphere median seeing (3free is 0.51", its mode is 0.35". The best 
seeing (minimal OT strength) is observed in October - November. The typical median value for that period is 
« 0.83". 

Median corrected atmospheric time constant is 2.58 ms, it exceeds 3.3 ms in conditions of weak turbulence. 
Median isoplanatic angle is 2.07" over whole period and 2.38" under condition of the seeing better than its 
median. Features of vertical distribution of the OT are presented in paper.'^ 




Figure 3. Layer-by-layer cumulative distributions of the OT above Mt. Maidanak. Left: processed with Turbina restoration 
algorithm Right: with help of the new algorithm 

4. REVISION OF THE OT PROFILE RESTORATION 

Turbina restoration algorithm is based on calculation of mean scintillation indices over accumtime (1 minute) 
and direct minimization of non-linear system in unknown J^^^. The obtained cumulative distributions doesn't 
look physically and may be explained with joint effect of the restoration errors and non-negativity restrictions. 
This effect leads to underestimated characteristic points (medians and low quartiles) of OT values in separate 
layers.'^ 

New restoration algorithm {atmos-2.96.0) is based on Non-negative least squares (NNLS) method developed 
40 years ago for the solution linear equations in terms of least squares with non- negative restrictions. New 
algorithm processes 1 s scintillation indices and then averages 1 s solutions. As result we have more accurate 
profiles and can use more layers (more than 12) in the restoration. The cumulative distributions of the OT 
intensity in 6 layers after reprocessing with atmos-2.96.0 are shown in Fig. [3] on right. 

Further modification of the profile restoration algorithm {atmos-2.96.7) was done in 2010 for the processing 
of 2-years set of measurements at Mt. Shatdzhatmaz. This version combines both DIMM and MASS data 
processing. Reasons are: 

• To calculate propagation effect in differential motion. 

• To restrict non-physical negative OT in ground layer. 

• To determine at once OT in ground layer. 

An example of such processing for 2009 February 19 is shown in Fig. S) For such modification we had to 
introduce DIMM weighting functions. 

4.1 DIMM weighting functions 

DIMM standard theory is based on near-field approximation what leads to partial loss of high turbulence power. 
To have an equations similar to MASS equations, the weighting functions (WF) were introduced for the differ- 
ential motion. The WFs provide next relation for transversal and longitudinal motion: 

= / clih)Wi4h)dh. (1) 
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Figure 4. Evolution of the OT above Mt. Shatdzatmaz in 2009 February 19 night as processed by atmos-2. 96. 7 version 



To calculate WFs we start from earlier studies of Martin^ and Tokovinin.l^ For example, WF of the transversal 
differential motion for the case of G-tilt can be written as 
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Polychromatic WFs can be computed by replacing the usual Fresnel term cos^(7rAft,/^) with the polychromatic 
Fresnel filter which is a square of the real part of Fourier transform of the incoming radiation energy distribution 
similarly it doing in MASS WFs calculation.'^ 

In Fig. [5] (on right) some WFs calculatcc0 for different DIMM devices are shown. Using the WFs technique 
very compact DIMM device is possible to be built. Of course, parallel measurements with MASS are needed to 
provide OT vertical distribution. 

4.2 DIMM finite exposure correction 

DIMM output files contain the correlation coefficient p between adjacent measurements of the image positions. 
The correction is bases on its calculated dependencies on p value. The method is approved theoretically. 

For 4 ms exposure and 200 frames/s, linear approximation of the correction is in the form: 



CT^ = ct2(i + o.15(p-1)). 



(3) 



* http: / / curl.sai.msu.ru / mass / download/doc / dimm_specs.pdf 
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Figure 5. Left: WFs for our DIMM (D = 0.09 m, B = 0.196 m), lAC DIMMA (D = 0.05 m, B = 0.20 m) and hypothetical 
mini-DIMM (D = 0.05 m, B = 0.05 m). Right: SAI DIMM WFs for Z-tilt (solid), G-tilt (dashed) and SR-tih'' (thin 
line). Longitudinal motion - black, transversal - red 



where is measured variance and cr^ — corrected 

The median of measured p = 0.85 what produces the correction w2%. Only 10% of data requires correction 
greater 5%. After December 2009 the exposure was reduced to 2.5 ms, so needed correction became twice less. 

5. WIND IN MASS DATA 

In MASS we see temporal fluctuation of the measured radiation mainly due to turbulence translation by the 
wind. In general, the variance of the flux fluctuation for the case of finite exposure r, 

<jI = J C^{h)W'{h,T,w)dh (4) 

where is the modified weighting function depending not only altitude h but wind profile w{h) and exposure. 

poo 
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where Ag is wind shear spectral filter having analytic expression.^i^ 

Two ultimate regimes can be emphasized: short exposure: wt <C D and long exposure: wt ^ D. In these 
cases we can take wind outside the weighting function integral. There are few problems related to this topic: 

• MASS finite exposure correction 

• Atmospheric time constant estimation 

• Potential photometric accuracy evaluation 

• Wind vertical profile extraction 

5.1 Potential photometric accuracy evaluation 

In the long exposure regime (t > 0.1 s for MASS apertures) 
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Figure 6. Left: Evolution of the median values of S3 index over entire period of Maidanak campaign. Right: Comparison 
of the wind w (filled dots) and wind at 70 mbar V70 (empty circles) 



From MASS data we can evaluate directly the index 6*3 introduced in paper Kenyon et al.'^ The index defines 
a photometric accuracy (scintillation noise) in aperture D averaged over exposure r: 



^2 = D-^/\-'Sl 

In the Figinion left, the evolution of the S3 index over 2 years of the Maidanak campaign is shown. 

5.2 High-altitude wind from MASS data 

In the long exposure regime some effective wind cam be calculated: 
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where averaging is done with weight C^{h)h^ which is maximal at the height of « 16 km above the summit 
corresponding « 70 mBar. In the FiglHlon right, the comparison of such value with data from NCEP/NCAR 
database is shown. 

The studies give us an assurance that on the base of extended MASS data (indices for set of exposures, which 
are collected last 2 years) the wind profile extraction is possible. Initial approximation may be easy obtained 
from long exposure indices, then non-linear equation set must be solved. 



6. SOME PROBLEMS 

6.1 The question No 1 

What does really DIMM measure: Z-tilt or G-tilt? The cost of this uncertainty reaches 12 - 17% in OT strength 
(see Fig. [5] on right). When these data are used for ground layer (GL) turbulence estimation, inaccuracy as large 
as « 30% can appear. 

Calculations show that erf /af = 1.61 for G-tilt and 1.51 for Z-tilt. Real measurements give the ratio erf /uf = 
1.56. High altitude turbulence can increase these ratios for both tilts. So z-tilt model is slightly preferable. 
Distribution of the measured erf /erf is quite wide, but corresponds an accuracy of the measurements, so data 
filtration by this ratio is never possible. 

We need reliable experimental verification method to check any instrument and any image centering algorithm. 




Figure 7. Left: Dependence of slow OT fraction on ground wind (circles), dashed lines are 1-st and 3-rd quartiles. Solid 
lines show calculated dependence in the case of Kolmogorov model. Longitudinal case is presented in black, transversal 
— in red color. Right: Ground layer turbulence versus OT in free atmosphere for two cases: full differential motion power 
(blue points) and high frequency differential motion (black points). 

6.2 The question No 2 

Our DIMM control software (dimm) produces the differential motion variances each basetime (1 - 2 s). Therefore, 
we can analyse the motions in two timescales: faster and slower than 1-2 s. We have detected that low frequency 
motions have a power much greater than Kolmogorov model predicts. For low motions the median erf / (Jt = 1-16. 

Hence, we observe the non-kolmogorov, low altitude, slow turbulence. Dependence on the ground wind is 
evident that is show in Fig [7] on left. Is it local turbulence which isn't connected with atmosphere itself? 

If we refuse this power then overall seeing decreases from 0.93" to 0.86". More important that GL turbulence 
diminishes in 1.3 times (see Fig [7] on right). Should we include low frequency power in DIMM results? 

6.3 The question No 3 

In the last version atmos-2.96.9 the strong scintillation correction was excluded to have the strong scintillation 
effect undistorted. Only exact conversion of Rytov variances to scintillation indices was used. It is not enough 
what one can see in Fig. [7] on right. The impact of strong scintillation is seen as sharp drop after 10^^^ m^^^ 
free atmosphere OT intensity. 

The correct estimation of this effect is more important in the range 10"^'^ . . . 10^^^ m^/'^, where it is unevident. 
Unfortunately, there is no clear theoretical description of the effect for real astronomical conditions. Probably, 
an iterative processing should be needed. 

6.4 The question No 4 

Atmospheric coherence time tq derived from MASS data is biased. Some ways are possible to correct it: experi- 
mental calibration,^^ wind profile integration, modification of existing DESI method.!^ 

It was presented already that we plan to develop the restoration of the wind profiles together with the OT 
profiles. As long as an accuracy of these wind profiles is unknown, we can not estimate resulting tq errors. 
Probably such way will be impractical. 

It may be shown that for short exposure regime wt <^ D\rF the next approximation is valid: 



al = al-^^ j Cl{h)w{hfU{h)<lh 



(9) 



The second term is the correction to zero exposure for certain index. The modification of DESI UK^thod means 
an usage of the full set of temporal covariances to obtain needed combination for atmospheric winds moment 
'W2- Rejection of used empirical calibration is assumed as well. 

7. CONCLUSION 

The main result of our researches as in the practical use of the MASS/DIMM device and well in the processing 
and analysis of data of the measurement is that there is not yet the potential of both methods fully realized. 
Clearly, a simple mechanical combination of these two devices is not enough, also data processing should be 
common. 

The posed above questions reflect the problems which need to be solved in the nearest future to provide more 
accurate and more robust data on optical turbulence in the earth atmosphere. 
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